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We study the electronic structure of a spherical jellium in the presence of a central Gaussian impurity. We test 
how well the resulting inhomogeneity effects beyond spherical jellium are reproduced by several approximations 
of density functional theory (DFT). Four rungs of Perdew's ladder of DPT functional, namely local density 
approximation (LDA), generalized gradient approximation (GGA), meta-GGA and orbital-dependent hybrid 
functionals are compared against our quantum Monte Carlo (QMC) benchmarks. We identify several distinct 
transitions in the ground state of the system as the electronic occupation changes between delocalized and 
localized states. We examine the parameter space of realistic densities (1 < rs < 5) and moderate depths of the 
Gaussian impurity {Z < 7). The selected 18 electron system (with closed-shell ground state) presents Id 2s 
transitions while the 30 electron system (with open-shell ground state) exhibits If 2p transitions. For the 
former system, the accuracy for the transitions is clearly improving with increasing sophistication of functionals 
with meta-GGA and hybrid functionals having only small deviations from QMC. However, for the latter system, 
we find much larger differences for the underlying transitions between our pool of DFT functionals and QMC. 
We attribute this failure to treatment of the exact exchange within these functionals. Additionally, we amplify 
the inhomogeneity effects by creating the system with spherical shell which leads to even larger errors in DFT 
approximations. 



PACS numbers: 

I. INTRODUCTION 

Kohn-Sham density functional theory''^ (DFT) is now ar- 
guably the most popular computational method in theoretical 
condensed matter physics and materials science. While the 
method is exact in principle, in practice it is applied only with 
approximations to the unknown, exact exchange-correlation 
functional. The conventional semilocal approximations of- 
ten fail to describe the structural, defect, and other properties 
of materials with strong electron-electron correlations — a fail- 
ure that is not only quantitative but often qualitative. Indeed, 
even for less strongly correlated materials where the method 
yields reasonable predictions, further increases in accuracy 
are highly desired. 

On the other hand, the quantum Monte Carlo (QMC) 
method allows direct solution of the many-body problem of 
interacting electrons using stochastic techniques^-S. In fact, 
the local density approximation is based on QMC calculations 
of the homogeneous electron gas^ (HEG). The only significant 
source of systematic errors in the QMC method is the fixed- 
node approximation^ to the Fermion sign-problem. Using 
fixed-nodes, QMC has proven to been very effective in pro- 
viding high accuracy results for many real systems such as 
molecules, clusters and solids with hundreds of valence elec- 
trons that are within 1-3% of experimen t^-'"'" . More recently, 
new techniques have been developed to reduce the fixed-node 
errorsi^^, further improving the accuracy of the technique. 

The spherical jellium system has been extensively in- 
vestigated as a model of large clusters of simple met- 
als (see, e.g., review ITsIl ). It has been found to have 
pronounced shell structure with magic numbers N = 



2,8,18,20,34,40,58,92,... approximately corresponding 
to fully-filled closed shell of each orbital momentum as 
ls^\lp^\ld^^\2s^\lf^^\2p^'\lg^^\2d^Hh^^3s\ A few bench- 
mark QMC studies existi^— for the total and correlation en- 
ergies of closed-shells at magic numbers with up to 106 elec- 
trons. These results have been later compared against sev- 
eral DFT approximation s ' . Surface correlation energiesi^ 
and surface exchange-correlation energie a'^-^" have been ob- 
tained using extrapolated values to infinite size. In general, 
the values obtained for correlation surface energies obtained 
with spheres are consistent with the latest DMC results for jel- 
lium slabs^ and other methods, but differ to some extent from 
the original DMC results from slab geometries^^. 

In this paper, we benchmark four levels of DFT approxi- 
mations, following the "Jacob's ladder" metaphor of Perdew 
and co-workers21"— , by performing accurate quantum Monte 
Carlo calculations on the model of the interacting electron gas 
subject to a central impurity potential. The interacting elec- 
tron gas is interpreted as a finite jellium sphere with N elec- 
trons of average density r^. An attractive spherical Gaussian 
of tunable strength at the origin represents the impurity poten- 
tial. The proposed system closely relates to an atom in a real 
material while retaining a simplicity that is amenable to highly 
accurate solution under a wide range of conditions (slowly- 
and rapidly-varying density regions as well as a wide range of 
density values). The purpose of solving this simple model is 
to understand the role of electronic correlations and exchange, 
to understand the features that more sophisticated functional 
must possess in order to describe, e.g., strongly coiTelated ma- 
terials, and to provide essentially exact benchmarks for test- 
ing new DFT approximations. Tabulated data of our results 
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is given in Electronic Physics Auxiliary Publication Service 
(EPAPS) Document No. []. 



II. MODELS WITH SPHERICAL SYMMETRY 

A. Spherical Jellium Model 

One of the simplest ways to neutralize the negative charge 
of N electrons is to consider a sphere of positive charge ps of 
uniform density 

with rs as adjustable parameter corresponding to a Wigner- 
Seitz radius in solids and Rc as the sphere radius. This is a 
basic description of the spherical jellium model. The external 
potential due to the positive background is then 
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The Hamiltonian of the electronic system (in atomic units) has 
the form 
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where we have introduced the constant Coulomb self energy 
of the positive background as 
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The addition of Esei / ensures that in the thermodynamic limit 
the energy of jellium spheres approaches the energy of HEG. 



B. Spherical Jellium Model with a Gaussian Impurity 

In order to further enhance the inhomogeneity effects in 
the jellium spheres, we propose the addition of the attractive 
Gaussian shaped impurity at the origin. The external potential 
is then modified as 



K^,(r) = -Zcxp(-r70 + 
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where Z and a represent the depth and the width of the added 
Gaussian. The Gaussian form for the impurity is not meant to 
exactly describe core levels of a real atom, but to create addi- 
tional localization for valence electrons. Moreover, using this 
model, we can also avoid the locality approximation required 
to evaluate conventional pseudopotentials^. As a welcome 



consequence, the Gaussian potential also preserves the cusp- 
less property of single-particle orbitals at the origin. To sum- 
marize, by controlling the amount of localization at the im- 
purity, we can alter the shape and occupation order of single- 
particle states. We note that a similar potential has been used 
in studies of hetero-atomic clusters of Refs. l27ll28i 



C. Spherical Jellium Shell with a Gaussian Impurity 

Another possibility to further amplify the inhomogeneity 
effects in the jellium is to create a spherical jellium shell. This 
can be achieved by combining external potentials [Eq. of 
a larger system with N+M electrons and smaller system with 
M electrons to form 

VHcuir) = y:Lr\r) - T//i(r) - Zexp(-rV^'), (6) 

where integer variable M < N controls the geometry of a 
shell [inner radius R^^ — rsM^^"^ and outer radius = 
rs{M + N)^/^]. Note that in Eq. ^ the Gaussian impurity at 
the origin is also explicitly included. Above external potential 
closely resembles the potential of hollow clusters-^'^^". The 
difference in our model is the addition of Gaussian impurity 
which attracts electron density towards the origin. In turn, 
the enhanced inhomogeneity (i.e., charge separation between 
center and shell) provides even more severe test for the single- 
particle methods studied here. 



III. METHODS 

The goal of this paper is to compare the total energies 
and radial densities of equivalent states of the spherical jel- 
lium systems in Hartree-Fock (HE), DFT and quantum Monte 
Carlo Methods. Due to the spherical symmetry of the sys- 
tem and because the spin-orbit interaction is not considered, 
the eigenstates of H must also be eigenstates of the angu- 
lar momentum operators and AIl, spin operators and 
Ms- The advantage of the spherical symmetry is that rigorous 
upper-bound theorems apply for QMC and ground state DFT 
can be generalizec^ for the lowest energy state states of each 
value of L and S. Therefore, meaningful phase diagrams can 
be constructed with each method. 



A. Criteria used to select eigenstates of L and S 

For benchmark purposes, we only need to compare equiv- 
alent states with DFT and QMC. For comparison, we have 
selected states that are Ukely but not guaranteed to be the min- 
imum energy configurations. Since we are concerned with jel- 
lium densities of 1 < Ts < 5, our system can be characterized 
as weakly interacting. In this regime it is a reasonable to as- 
sume that, in analogy with atoms and 3D quantum dots-'^, the 
filling of the orbitals within a single shell of spherical jellium 
follows Hund's rules, i.e., for a given electron configuration, 
the state with maximum multiplicity (25" + 1) has the lowest 
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energy, and for a given multiplicity, the state with the largest 
value of L is likely to have the lowest energy. It is straightfor- 
ward to construct the eigenfunction of ^•^+^L symmetry (in 
the Russell-Saunders term symbol notation) as the Slater de- 
terminant of single-particle orbitals with occupancies chosen 
such that Ml = Yli rrii ^ L and Ms = Yli = ^- For 
Ml 7^ 0, the determinants will be complex-valued, an issue 
important in the QMC context and discussed further bellow. 

The original second Hund's rule only applies to 'iD 
spherically-symmetric systems, since it involves the total an- 
gular momentum L. There is debate, however, on the exis- 
tence and form of a second Hund's rule, with redefined shells, 
in the case of parabolic models of quantum dots in 2D at the 
highly correlated limit (see, e.g., Refs. [33lf34l) . While that de- 
bate is interesting, it is beyond the scope of this paper to vali- 
date the Hund's rules in the range of parameters explored with 
our 'iD models. 



B. Single-particle Methods 

To obtain the quantities of interest within single-particle 
methods we employ a non-relativistic atomic solver, orig- 
inally implemented in Ref. jssll and currently part of the 
QMCPACK suite^^*", modified to handle pure HF method, 
LDA, GGA, meta-GGA and hybrid functionals. The single- 
particle orbitals are conveniently expressed as radial-functions 
Rn.i (r) multiplied by an angular-part given as a spherical har- 
monic Y;, 4>). This choice ensures well defined quantum 
numbers. 

The radial equations [Rn,i (r)] are iteratively solved via Nu- 
merov algorithm on a large logarithmic grid (~ 5000 points) 
until the tight convergence criteria are met (AEtot < 10^^ 
and AEeig < 10~^^). For the later cases of the open shells 
we use the spin-unrestricted formulation of the HF and DFT 
theories. In short, in our formulation, we perform the spin- 
dependent average of the effective potentials within each sub- 
shell of the same n and /, resulting in a single radial function. 
As we discuss in the next subsection, these radial functions 
are then directly imported into QMC methods. 

In general, the Kohn-Sham and HF wave functions have 
a form of a single Slater determinant and are explicit eigen- 
states of Ms, but they are eigenstates of 5^ only if has 
the maximum value. Therefore, the L and S eigenstates with 
maximum multiplicity (25 + 1) and with the largest value of 
L automatically have the correct spin symmetry. 

In order to simplify the DFT implementation, for open 
shell cases, we have neglected the angular dependence of 
the spin-densities for calculation of the exchange-correlation. 
This treatment is commonly known as the spherical approx- 
imation used in many DFT atomic solvers for production of 
pseudopotentials. A very good numerical agreement (bet- 
ter than 1 mHa) was achieved between our and other atomic 
solvers (OPIUM^., FHWSPP^^'^ and APB^) for several open 
and closed shell atomic states and functionals. The four rungs 
of Perdew's ladder of approximate DFT functionals are rep- 
resented by the Perdew-Wang (PW) LDA^S, Perdew-Burke- 
Ernzerhof (PBE) GGAii, Tao-Perdew-Staroverov-Scuseria 



(TPSS) meta-GGA24 and PBE hybrid (PBEO)^ exchange- 
correlation functionals as implemented in the LIBXC library^ 
and in the Quantum Espresso suited. 

C. Quantum Monte Carlo Methods 

The trial many-body wave function serves as the most im- 
portant input for the quantum Monte Carlo methods within 
the fixed-node diffusion Monte Carlo approach. For this prob- 
lem, following the approach used in previous studiesi^^, we 
choose a many -body wavefunction that is a product of spin-up 
and spin-down Slater determinants (D) and a Jastrow correla- 
tion factor (J), written as 

*T(ri,r2, . . . ,r7v) ^D[ipl{ri), . . .,ipl.j{rM)] 

X D[ip\{rM+l), . . . , fN^Mi^N)] 

X exp[J(ri,r2, . . . ,rAr)], (7) 

assuming the first M electrons to be spin-up and the remain- 
ing — M to be spin-down. The Slater determinants are 
constructed from orbitals generated in single-particle methods 
and have a form 

The symmetric Jastrow correlation factor includes well- 
known electron-electron cusp conditions as well as one 
and two-body correlation functions (for details see, e.g., 
Ref. lHH). As a note, we did not find it necessary to include 
the multipolar terms into the Jastrow factor as in Ref. [l^ 

The Jastrow term is further variationally optimize d'^'^^ . 
The optimal 4't is then used in diffusion Monte Carlo (DMC) 
method to obtain the ground state fixed-node (FN) energies 
and other expectation values of the system. As a additional 
step, for selected states, we also perform reptation Monte 
Carlo (RMC) calculationsi^i to obtain pure expectation values 
of the fixed-node densities. All the QMC calculations were 
performed using QWalk code^. 

As we have already mentioned, the L and S eigenfunctions 
considered in this paper with Ml 7^ are complex-valued and 
require the use of the fixed-phase DMC algorithm^. How- 
ever, the associated fixed-phase errors are in general different 
(and possibly larger) than the fixed-node errors for real- valued 
wavefunctions^-. To avoid mixing results with two different 
approximations, we exclusively use the real-valued Ml = 
projections of the same LS eigenfunction (as they are ener- 
getically degenerate). The Ml = projections can be readily 
obtained from Ml = L states by the recursive application of 
the momentum lowering L~ operator. As a consequence, the 
complex-valued single Slater determinant is replaced by the 
real-valued linear combination of Slater determinants. Please 
refer to the Appendix for the linear combinations used for 
each open shell state. 

IV. RESULTS 
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A. Tests and validations 

To test the numerical implementation of our HF, DFT and 
fixed-node (FN) DMC methods we have recalculated previ- 
ous results for the closed-shell jellium spheres. We achieve 
an excellent agreement with the PW LDA and PBE GGA 
ener gies of Ref. [l9| and with TPSS meta-GGA energies of 
Ref . uOi Obtained results for HF energies are identical to 
Ref.|5l]for model sodium clusters of rg =4.0 with up to 196 
electrons. However, we find lower HF energies (by 0.5 mHa 
per electron) when compared to energies of Ref. J_8s. This is 
most likely due to the non-self consistent treatment of HF in 
Ref. in Finally, all our FN-DMC energies using LDA orbitals 
are within error-bars of FN-DMC energies of the Ref. [HI The 
sole exception is the 2 electron system at rg = 1, where the 
differences are larger. 

The single limitation for the DMC calculations comes from 
the use of fixed-node approximation. Introduced fixed-node 
errors can be reduced by expanding the trial wave function 
in multi-determinantsi^iii or by using the backflow correla- 
tion correctionsSi"— , or both. Previously, it was found that 
the QMC calculations for the closed shell jellium spheres'^ 
did not suffer from large fixed-node errors (by comparing to 
FN-DMC results with small multi-determinant expansions'^). 
In order to systematically check for these errors, we have ex- 
tended the previous FN-DMC calculations to include back- 
flow correlations for all the densities and particle numbers. 
We find small and uniform gains to correlation energies on 
the order of 0.3 mHa per electron. Therefore, it is very rea- 
sonable to assume that, for the energy differences considered 
here, these small corrections cancel out. 

B. Results for the spherical jellium with impurity 

1. Id — >■ 2s transitions for N= 18 

The excellent agreement for the closed-shell energies of jel- 
lium spheres, obtained with both with QMC and DFT meth- 
ods, with earlier publications gives us confidence to study 
open shells. As the first testing case, we have selected a jel- 
lium sphere with 18 electrons, which for the impurity free sys- 
tem is one of the fully-filled closed shells and corresponds to a 
cluster with a stable configuration. As we increase the attrac- 
tive Gaussian potential of the impurity, the nearest unoccupied 
2s level is lowered in energy below the Id level. Interest- 
ingly, exactly the same crossing of the levels was previously 
observed in the context of hetero-atomic clusters2i2^ and 3D 
quantum dots'*-. As a result, the ground state of the system 
will change from closed shell occupation ^ S{ls'^lp^ld^'-*) to 
either ^D{ls^lp'^ld^2s^) or ^F{ls^lp^ld^2s^) open shell 
occupations. As a note, each state differs in occupation in 
more than one shell {Id and 2s) while Hund's rules strictly 
apply only to occupations within a single shell. 

The phase diagram depends on the variables r^, Z and a. 
We limit our discussion to realistic density range between 
rg = 1 and rg = 5 found in the majority of bulk systems. 
Since the <t and Z parameters are effectively coupled together. 



TABLE I: Comparison between the total energies of the closed-shell 
and open-shell states for TV = 18 and = 1. The latter state 
is clearly lower in energy. 

1^ E^MC ^LDA ^PBE 

0.39104(2) 0.39317 0.38908 0.42786 
0.38817(2) 0.39030 0.38628 0.42386 

we choose to fix the width of the Gaussian at cr = 0.6 and vary 
only its depth Z. As a result, we find robust Id 2s transi- 
tions within < Z < 7. We summarize our calculations in 
the phase diagram, Fig.[T] 

There are several general observations which can be drawn 
from the phase diagram, Fig.[T] As we increase the attractive- 
ness of the impurity, we see changes on the occupation of a 
more delocalized Id orbital to more localized atomic-like 2s 
orbital in the following order: Id^" ld^2s^ ld^2s^. 
In addition, we notice the existence of high density region 
(vs < 1.4) which has always partial 2s occupation (see also 
Table This region for jellium spheres at — 1 was not 
discussed in previous studies. 

From a methodological point of view. Fig. [T]provides de- 
tailed comparison between the single-particle methods and 
our DMC benchmarks. Our first observation is that HF greatly 
overestimates while LDA underestimates the region of the sta- 
bility of the ^D{ls^lp^ld^2s^) state. In fact LDA and HF 
bracket our best estimate, which is not surprising considering 
the success of hybrid functionals. Second, the GGA rung of 
the functional ladder represented by the PBE leads to clear 
improvement in accuracy over LDA. Next, the TPSS meta- 
GGA and PBEO hybrid functionals agree even closer with 
DMC predictions than PBE however without particular order 
in accuracy. 

Last, to give some measure of the fixed-node errors, we also 
compare the DMC results using LDA orbitals with DMC re- 
sults using HF orbitals. Despite the different nodes resulting 
from these choices the DMC energies are very similar indicat- 
ing that the nodal eiTors remain small. 



2. I f — !> 2p transitions for N=30 

In general, the behavior of the 30 electron system with a 
partially-occupied 1/ shell is similar to the case with 18 elec- 
trons. The main difference is that, for the impurity free sys- 
tem, the ground state configuration is characterized by a high- 
spin quintet state at I f shell and with 2p as the closest empty 
level. Therefore, the role of the extended orbital that transfers 
the electron is taken by 1/ while Ip localizes and receives an 
electron as the impurity potential is applied. 

Since Ip is more degenerate than Is, the 30 elec- 
tron systems allows us to study more complicated 
effects of correlations and exchange at the impurity 
as we change the impurity potential. The accessi- 
ble states of the interest are ^/(ls^lp*'ld^°2s2l 
7/(ls2l/ldi02s2l/92pi), 9G(ls2l/ldi02s2l/82p2j and 

"^(Is2lp6ldl02s2ij72p3)^ 
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FIG. 1: Phase diagram of spherical jellium with 18 electrons and 
for fixed a — 0.6 as a function of density Vs and impurity poten- 
tial strength Z. The transitions between the three lowest lying states 
are shown for unrestricted Hartree-Fock (HF) (red dot-dashed), PW 
LDA (blue dashed), PBE GGA (green short-dashed), TPSS meta- 
GGA (purple dot-dashed), PBEO hybrid (orange dot-dashed), DMC 
with UHF orbitals (black dot-dashed) and DMC with LDA orbitals 
(black full) lines (see also text). The lines in the figure are extrapola- 
tions over the discrete points by cubic splines. 

In Fig.|2]we directly compare the phase diagram obtained 
with DMC and mean field methods used in the 18 electron 
case. The detailed comparison shown in Fig. |2] reveals sev- 
eral important features. First, LDA, PBE GGA and TPSS 
meta-GGA, all semilocal functionals, predict almost identi- 
cal boundaries. When compared with our DMC results, we 
find that only the locations of ^/ <r^^^ S and o / phase 
transitions agree well, while the ^/ <->^ G and -f-)-^^ S 
transitions are shifted to lower Z values. Second, PBEO hy- 
brid corrects slightly for the lower Z shifts in the ^/ G 
and ^G S transitions. In contrast, HF method produces 
much more satisfactory agreement with DMC for smaller den- 
sities {rs > 2.5) but greatly overestimates the stability of ^^S 
state at higher densities (j-g < 2.2) due to the missing corre- 
lation. The behavior described above suggest that a) the inho- 
mogeneity effects in the density are relatively small and well 
captured by the semilocal functionals; b) the full non-local ex- 
change absent at the semilocal level and only partially present 
(25%) in PBEO functional is needed to correct for lower Z 
shifts. 



C. Results for the spherical jellium shell with impurity 

1. Id ^ 2s transitions for N= 18 

As in the first case of the spherical jellium with impurity, 
we choose to study the 18 electron system in the spherical 
shell potential [Eq. ^]. The occupation of the single-particle 
states for hollow cluster is assumed to be the same as for jel- 
lium spheres (as confirmed by Ref. [s^)- Rather than finding 
the extensive phase space diagram of the system we limit our 




1.5 2 2.5 3 3.5 4 

FIG. 2: Phase space diagram for jellium with 30 electrons and for 
fixed a = 1.5 as a function of density and impurity strength Z. 
We have identified ^I, and "5" states as occupying this sec- 

tion of the rs-Z space. The convention for lines is identical to Fig.[T] 



study to a very small subset of the space. We select M = 9, 
rs = 3.0 and Z ^ 8.5 to illustrate a size of localization errors 
from HF and DFT based theories. 

Figure [3] compares the energy difference between and 
states as a function of Z calculated with single-particle 
methods and DMC. Arguably, the main point of interest in 
Fig. [3] is the slope around F transition and secondly 

the size of the relative errors. As in previous subsection, all 
semilocal functionals (i.e., LDA, PBE, TPSS) results have a 
very similar slope (with half the steepness of the DMC curve). 
On the other hand, the slope of the HF curve is almost identi- 
cal to DMC curve. Not surprisingly, PBEO partially recovers 
the correct slope as it contains 25% of exact exchange. 

Finally, it is also instructive to analyze the radial densities 
for the and ^F states (see Fig. |4|i. As a benchmark we 
use the pure expectation of the density operator from reptation 
Monte Carloi2.(RMC) method with LDA orbitals. The densi- 
ties for each state and spin channel in Fig. |4] have a distinct 
double peak structure - the smallest is due to the presence of 
the Gaussian impurity and the largest due to the spherical shell 
itself. Also visible is the relative reduction of the smaller peak 
for the '^D spin-down channel due to the absence of the more 
localized 2s state. 

From the direct comparison with RMC results we deduce 
that the density at inner shell regions (i.e., smaller peak) is 
better described within HF while LDA and PBE GGA pro- 
vide better densities for the outer shell regions (i.e., larger 
peak). The PBEO hybrid smoothly interpolates between PBE 
and HF limits with the best overall description. Lastly, the 
TPSS meta-GGA results in an substantial increase of density 
at both peaks and decrease at the tails when compared to LDA 
and PBE GGA. We find this overcompensation of the den- 
sity, presumably due to the kinetic energy density contribu- 
tions specific to the TPSS functional, to be surprisingly high. 
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FIG. 3: Total energy difference between D and F states of the 
spherical jellium shell of 18 electrons with M = 9, ra = 3.0 and 
a = 0.6 in a small window of Z. All the methods are described in 
the text. 

V. SUMMARY 

In conclusion, we have studied the performance of a variety 
of DFT functionals with increasing complexity against quan- 
tum Monte Carlo benchmark in a spherical jellium model and 
an spherical shell model with an attractive Gaussian impurity 
at the center. The tunable strength of the impurity allowed 
us to find a number of interesting transitions between closed 
and open shell states. Our results shows that the development 
of better density functional approximation is increasingly re- 
quired as the system departs from the perfect spherical jellium 
case. We report several regions where the employed approxi- 
mations of DFT fail to find the same ground state as identified 
with QMC methods. 

The Id 2s transitions in 18 electron system are well 
described on the highest semilocal level (TPSS meta-GGA) 
as well in global hybrid (PBEO) DFT. On the other hand, not 
all 1/ — > 2p transitions in the 30 electron system were ac- 
curately captured. We argue that even as the inhomogeneity 
effects in the density are relatively small and well captured by 
the semilocal functionals, full non-local exchange is needed 
to accurately describe the system. Our work therefore further 
supports the need for the hyper-GGA functionals^ with fully 
non-local exchange and accompanying balanced correlation. 

In the spherical jellium-shell model with an impurity at the 
center, where the inhomogeneity in the electronic density is 
increased, the DFT methods with exact-exchange give better 
agreement for the studied transitions. The radial electron den- 
sities in the inner region closest to the impurity are correctly 
described at the HF level, while LDA and PBE GGA are more 
accurate in the outer region. The PBEO global-hybrid results 
smoothly interpolates between HF and GGA. We find surpris- 
ingly high deviations in density for the TPSS. Finally, we also 
publish our results in the EPAPS Document No. [] with the 
purpose of allowing a detailed comparison with newly devel- 
oped DFT functionals. 
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Appendix: Real-valued LS eigenfunctions 

Application of {L~)^ operator on the LS eigenfunction 
with Ml = L leads to real-valued LS eigenfunction with 
AIl = 0. The linear combinations of Slater determinants for 
the open-shell states in 18 electron system are then 

=Di (9) 

^_Lp2 + 2i^3), (10) 

where L>i = (2", 1", -1", -2"), D2 = (1^,0", -1") and 
D3 = (2^, 0^, —2^) are determinants with indicated rele- 
vant occupied orbitals (numbers stand for the orbital quantum 
numbers and superscripts indicate the spins). 

The Ml ~ linear combinations for the 30 electron system 
are 



*[«G(l/^2p2t)] 



1 

^7f 



V2{Di + D2) + 



(11) 



where = (~l-,l+,0+), = (l-,0+,-l+), D3 
(0-,l+,-l+) and 



n'L{if^2p^)] - 



1 



:[30i:'4 + 24155 -f6L»6 



6%/7\/TT 
+ 5V6D7 + 9\/5L»8 + sVSL'g 

+ 5V6D10 + 9V5Dn + 5VSD12], 

(12) 

where i:»4 = (1^,-1", 0+), D5 = (2^,-2-, 0+), Dq = 
(3-,-3-,0+), = {0-,-l-,l+),Ds = (l-,-2-,l+), 
D9 = (2-,-3-,l+), i?io = (l-,0-,-l+), Dn = 
(2-,-l-,-l+), = (3-, -2-, -1+) and 

1 



*r/(i/ 



V2V3V7VTI 

+ 5V2{Di6 + Dij 



[5Di3 + WDii + 9Di5 



(13) 



where D13 = (l-,0-,-l-), Du = (2-,0-,-2-), 
Di5 = (3-,0-,-3-), Die = (2", 1", -3-),i?i7 = 
(3^, —1^, —2^). Above results have been also verified nu- 
merically using code from Ref. is^- 
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r [a.u.] r [a.u.] 

FIG. 4: Radial densities for the spherical jellium shell of 18 electrons for (upper) "^F (lower) states with M — 9,rs = 3.0, Z = 8.5 and 
a = 0.6. Spin-up channel (left figure) and spin-down channel (right figure). All the methods are mentioned in the text and the convention for 
lines is identical to Fig.[T] The pure expectation of the density operator from RMC method employed LDA orbitals. 
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